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Supercurrent reversal in Josephson junctions based on bilayer graphene flakes 
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We investigate the Josephson effect in a bilayer graphene flake contacted by two monolayer sheet 
deposited by superconducting electrodes. It is found that when the electrodes are attached to the 
different layers of the bilayer, the Josephson current is in a tt state when the bilayer region is undoped 
and in the absence of vertical bias. Applying doping or bias to the junction reveals tt — 0 transitions 
which can be controlled by varying the temperature and the junction length. The supercurrent 
reversal here is very different from the ferromagnetic Josephson junctions where the spin degree 
of freedom plays the key role. We argue that the scattering processes accompanied by layer and 
sublattice index change give rise to the scattering phases which their effect varies with doping and 
the bias. Such scattering phases are responsible for the tt — 0 transitions. On the other hand if both 
of the electrodes are coupled to the same layer of the flake or the flake has AA stacking instead of 
common AB, the junction will be always in 0 state since layer or sublattice index is not changed. 

PACS numbers: 72.80.Vp,74.50.+r,74.45.+c,85.25.Cp 


I. INTRODUCTION 

Starting from a decade ago, two dimensional (2D) 
atomic layers synthesized and received a huge amount of 
interest.^ Graphene, the leading 2D material, has been 
studied a lot mostly because of promising applications 
in electronics, chemistry, optics, etc. beside unexpected 
relativistic-like electronic dispersion arising interest from 
fundamental point of view.^ It has been shown that 
although the low-energy quasi-particles in single layer 
graphene are massless Dirac fermions, however, the sit¬ 
uation for bilayer graphene (BG) is fundamentally dif¬ 
ferent revealing chiral gapless excitations with quadratic 
dispersion rather than linear.Subsequently electronic 
properties of mono- and bilayer graphene are far from 
each other and as the most famous example due to the 
so-called Klein tunneling the backscattering is absent in 
monolayer graphene (MG) while in the case of bilayer, 
the forward scattering is impossible.^ Very intriguingly 
one can simply use a perpendicular electric field to gen¬ 
erate a controllable gap in graphene bilayers.^^^^ 

Recently, it has been experimentally revealed that in 
epitaxial graphene, it is very probable to have step-like 
monolayer/bilayer (ML/BL) interfaces as well as bilayer 
flakes connected to monolayer regions.A natural 
question arises immediately about the interplay of mass¬ 
less and massive dispersions in the junctions containing 
both mono- and bilayer graphene regions. Earlier theo¬ 
retical investigations had proven that in the ML/BL in¬ 
terfaces the transmission probabilities can show a valley 
dependent asymmetry which suggests their usage in the 
generation of valley polarized electron beams.More¬ 
over in the presence of perpendicular magnetic fields the 
emergence of Landau levels with peculiar transport prop¬ 
erties has been studied. In particular, a rich Lan¬ 
dau spectrum has been predicted^^ and an asymmetry 
in the dependence of transport features to the sign of 
magnetic field and charge carriers has been experimen¬ 
tally observed.Other theory works have focused on 
the edge states properties and quantum transport via 
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Figure 1. (Color online) Josephson couplings: supercurrent 
streams through the biased BG flakes embedded between two 
Josephson coupled superconducting MG leads, Sl and Sr. 
The length and width of the junction is indicated by d and 
W, respectively. A voltage difference Vb is applied perpendic¬ 
ularly between the BG layers. The superconductor either are 
attached to top and bottom layer of the BG flake (a) or they 
are connected to the same layer of the BG flake (b). Qualita¬ 
tive description in top of each setups demonstrates the trans¬ 
mission of the quasiparticles through the BG flakes by the red 
and blue curved arrows. 


channels introduced by the interface.On the other 
hand, theoretical investigations of the transport through 
bilayer flake sandwiched between two monolayer nanorib¬ 
bons have revealed that the conductance oscillates be¬ 
tween maximum and zero as a function of bilayer flake 
length. 

Pioneering work of Beenakker^^’^^ and successive ex- 
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perimental realization of superconducting proximity ef¬ 
fect and supercurrent passing through graphene^^ proved 
that graphene-based superconducting heterostructures 
can have very interesting properties^"^”^^ which some 
of them are experimentally realized.Among the 
variety of theoretical predictions, the so-called specu¬ 
lar Andreev reflection, despite several proposals, still 
waits to be detected.In addition, when the nor¬ 
mal graphene region is replaced with a magnetic one, 
the well known 0 — tt transitions have been predicted to 
take place with some features different from conventional 
superconductor-ferromagnet-superconductor (SFS) junc¬ 
tions. Intriguingly it has been suggested that one can 
see the transitions by varying the doping of the magnetic 
region sandwiched between superconductors. 

In this paper we investigate the Josephson effect 
through a BG flake embedded between two monolayer 
sheets as can be seen in Fig. 1. We find that the first 
setup in which the two ML leads are coupled to the lower 
and upper layers of the flake, is in the Josephson tt state 
when the doping /i and vertical bias Vq are zero. Apply¬ 
ing gate voltages which results in doping and bias to the 
flake the device can undergo a tt to 0 transitions. More¬ 
over at the presence of finite /i or Vb, the 0 — tt transitions 
occur with temperature and the junction length as well. 
It must be noted that the reason of such transitions are in 
contrast with those taking place in the Josephson devices 
with ferromagnetic weak links. Here the origin of tt state 
is the phase factors of transmission coefficients through 
the flake which influence the superconducting phase de¬ 
pendence of the current carrying Andreev bound states. 
The appearance of phases is somehow related to the fact 
that quasiparticles passing through the flake in the setup 
(a) need to change their layer and sublattice index while 
in the setup (b) for all the modes the scattering phases 
are zero. Varying the bias voltage and doping of the flake 
changes the scattering properties of the flake and at cer¬ 
tain points the overall effect of the scattering phases are 
diminished leading to the 0 state Josephson effect. 

This article is organized as follows. In Sec. II, we intro¬ 
duce our model and explain the method we use to calcu¬ 
late the supercurrent that streams between two Joseph¬ 
son coupled superconducting leads. In Sec. Ill, we present 
and discuss our numerical results that shows the effect 
due to bilayer step junction as well as bias voltage on 
supercurrent. Finally, in Sec. IV, we conclude and sum¬ 
marize our main findings. 


II. MODEL AND BASIC FORMALISM 

We investigate the supercurrent flowing between two 
superconducting reservoirs on top of MG regions which 
are connected through a weak link containing a biased 
BG flake. In fact, there are two different configurations 
for the setup depending on how the BG flake is connected 
to the MG regions in the left and right. The left and 
right reservoirs can be attached to the top and bottom 


layers of the bilayer, respectively as Fig. 1(a) or both 
of the monolayer reservoirs are connected to the same 
layer of the flake as Fig. 1(b). To be precise, we will as¬ 
sume the the MG regions are fully covered by supercon¬ 
ductors and only the middle BG flake is in the normal 
state. We further consider an s-wave order parameter, 
A(r) = in the left and right superconducting 

regions and A(r) = 0 inside the flake, where Aq is an 
isotropic gap. 

Similar to any other weak link, the supercurrent across 
the junction is mainly carried by the discrete bound 
states which are the result of closed loop motion of the 
quasiparticles with subgap energies between two super¬ 
conductors, known as the Andreev bound states (ABS). 
One can see that it is usually sufficient to find ABS en¬ 
ergies in order to calculate the supercurrent carried by 
them. In the short junction limit, where the length of the 
junction (d) is much smaller than the coherence length 
of the superconductor (^), the current establishes due to 
the phase dependence of ABS energies and the continuous 
spectrum does not contribute since its density of states 
is almost phase independent. The relation between the 
Josephson current I passing through the junction with 
the transverse width W and subgap quasiparticles (e^) 
at finite temperature T is given by, 

r 2eVF /* „ V- . 

1 = - — / dky 2 _^ tanh 


2kBT \ dcj) ^ ^ 


The sum is over all ABS energies which are positive cor¬ 
responding to the different transverse momenta ky. 

The ABS energies can be calculated in the framework 
of Bogoliubov-de Gennes (BdG) equations which describe 
the superconducting correlations between particles and 
their time-reversed counterparts (holes). Exploiting the 
BdG equations, the ABS energies will be obtained as the 
roots of a characteristic equation containing the whole 
scattering matrix of the junction. This method is based 
on the fact that the ABS can be assumed as the states 
scattered completely to themselves inside the junction. 
The whole scattering matrix consists of two parts: S^ 
which describes the normal scatterings of electrons and 
holes within the bilayer flake and S^ responsible for the 
scattering away from the normal-superconducting (NS) 
interfaces. Since for subgap energies 5 < Aq there are 
no propagating modes in the superconducting regions L 
and R then the scattering matrix S^ can be defined with 
relation as = Sa^b which relates the outgoing to incom¬ 
ing quasiparticles. This matrix consists of four blocks in 
electron-hole (Nambu) space corresponding to different 
processes as, 

. ( 2 ) 

\^he ^hh J 


In general each block of the scattering matrix in Nambu 
space denoted by S have the following form. 


f ^ll tLR\ 

\tRL Tbr) ’ 


( 3 ) 
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consisting of reflection r and transmission t matrices 
which here are 2x2 matrices in the space of two lay¬ 
ers. The scattering matrix of the NS interfaces involves 
only the reflection processes and subsequently the corre¬ 
sponding normal and Andreev reflection parts are block- 
diagonal as below, 


S 


ee — 




(4) 


Now we should remind that at the left and right NS in¬ 
terfaces, depending on the configuration, either different 
layers (setup a in Fig. 1), or only one of the layers (setup 
b in Fig. 1) are involved in the reflection processes. As 
a result the matrices for the two setups are obtained re¬ 
spectively as below. 


(a): 

^ee 

= -a(e)f„ 

> = /3(e)e*‘^/2f„, 



Sfe 

= a{e)Td , 


(5) 

(6): 

Se^e 

= -a{s)fu 

> = /3(e)e*‘^/2f„, 



qB 

^ee 

= a{e)Tu , 

= /3(£)e-*^/2f„, 

(6) 


with Tu^v = ('^o A in which being the 

Pauli matrices in the layer space and ^hh = —See, 
SL = exp(-i(/))Sf';„ = exp(i(/))Sf;,. The normal 

and Andreev reflection amplitudes are obtained for 
the MG-based NS interfaces in Ref. 31 as 0 ^( 5 ) = 
iCsin^/[(5/Ao) cos^+C] and p{e) = cos6>/[(5/Ao) cos6>-l- 
C] in which C = yj (^/Aq)^ — 1 and 0 is the incidence an¬ 
gle. 

The scattering matrix of the normal BG region S^ re¬ 
lates two coefficient vectors of transmitted and reflected 
as Bb = Sbclb- Owing to the fact that matrix Sb does 
not couple the electrons and holes together, it has a block- 
diagonal form and given by. 


Sb 


/S(e) 0 \ 

V 0 ’ 


(7) 


Here S(5) and S*(—e) are the unitary and symmetric scat¬ 
tering matrices governing the scattering properties of the 
electrons and holes. The reflection and transmission ma¬ 
trices f and t for the BG flake embedded between two 
MG are given by. 


(а) : 

tll = riiTu 
^ / 
thR = G 2 ^+ 

( б ) : 


rRR = r22Td 
tRL = ^ 21 'f- 


tll = riiTu , 
tLR=tiiTu , 


rRR = r^^Tu 
tRL = tiifu 


(8) 

(9) 


that in the setup (b) the layer index is conserved, however, 
in the setup (a) since the only way to pass through the 
scattering region is via the interlayer hopping between 
the two layers it is not conserved anymore. Transmis¬ 
sion tij and reflection amplitudes can be obtained by 
matching the wavefunctions at the ML/BL boundaries. 
In order to complete the construction of the layer resolved 
scattering matrix, the excitation spectrum of BG as well 
as MG and their wavefunctions as the scattering basis 
are needed. 

The low-energy Hamiltonian of BG flake in the pres¬ 
ence of layer asymmetry due to the bias voltage applied 
perpendicularly and in the vicinity of non-equivalent cor¬ 
ners of Brillouin zone K and K' is given by. 


'Hbg 


/-/i + ^ vfP- t±_ ^ \ 

vfPf + y 0 0 

t± 0 vfPf 

\ 0 0 VfP- - y / 


( 10 ) 


in the layer and sub-lattice spaces with eigenfunctions 
of the form = (c^pp, c^pp, cb 2 ,p, ca 2 ,p). Here p± = 
Px ± ipy and p = (px^Py) is the two dimensional momen¬ 
tum measured relative to the K point, vf ~ lO^m /5 is 
the Fermi velocity and fi being the chemical potential. 
The nearest neighbor atoms in two layers A1 and B2 
are connected by interlayer hopping term t± 0.3 ev 
which tends to equalize the charge densities in the two 
layers. The potential difference between the two layers 
is involved by the parameter V which opens a gap in 
the spectrum in contrast to the case of gapless spectrum 
V = 0. The vertical bias through the BG also works 
against the interlayer hopping since it generates charge 
imbalance between the two layers. Taking plane wave ba¬ 
sis we end up for the excitation spectrum of the Hamil¬ 
tonian with four energy bands given by. 


(eBG + = (VFP? + (VV4) + {tl/2) 

±^{tl/2r + {vFPnV^+tl), ( 11 ) 


the corresponding eigenvector reads, 

= .4 (7- VfP+ V ^vppFj , (12) 

where r] = [ 7 ?. — {vFp)‘^]/t± with = sbg + 

fi ± V/2 and the normalization factor is M = 

[ivFPx (7- - (^BG + 

In order to finds the scattering matrix Sb we need 
the eigenfunctions inside the left and right MG regions 
which are immediately coupled to the superconductors. 
The wave-functions are represented in the space of two 
sublattices as = (cy^iy, CBi,k), and the two dimen¬ 
sional Dirac Hamiltonian around K and K' points reads. 


^MG = 


fi FF{kx tky^ 

VF{kx-\-iky) fi' 


(13) 


for the two different setups, respectively, with t± = 
(tx ± iTy)l2. The labels 1 and 2 denote the two layers 
of the BG flake. Therefore it is clear from above relations 


Here the two-dimensional momentum is'k = {kx^ky) and 
fi' is the chemical potential in the MG so that the exci¬ 
tation energy given by smg = fi' ^VFk is measured with 
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respect to the Fermi energy in MG. The corresponding 
eigenvector are, 

= ( 2 cos 6 »)”^/^ (1 ( 14 ) 

with ^ = L^R and 0 = diTctdin{ky/kx). Since the eigen¬ 
vectors will be utilized as scattering basis they need to 
be normalized in a way that each state carries the same 
amount of quasiparticles current density. 

We are now in a position to write down the wavefunc- 
tions in the left, middle and right region as a linear su¬ 
perposition of constructed scattering basis Eqs. (12) and 

(14), 


^kxX-\-kyy) 

_^^^i-PixX-\-kyy) ^^i{-P2xX-\-kyy) 

'^R = ( 15 ) 


here + (—) denote right (left) moving quasiparticles with 
r and t being the reflection and transmission amplitudes 
respectively. Scattering coefficients can be computed by 
solving the linear system constructed from the matching 
boundary conditions at the interfaces and subsequently 
will be used to find the the transport properties across 
the junction. Thus we need to provide the appropriate 
boundary conditions at the boundaries x = 0 and x = d. 
The MGs can be attached either to the same layer of BG 
or to the different layers. We assume the zig zag bound¬ 
ary conditions at the ML/BL interfaces. The boundary 
conditions at the left interface x = 0 are. 


«'l(x = 0 )U, =«'b(x = 0 )|^^ 

«'L(a: = 0)|Bi =^'B(a: = 0 )|B^ 

«' B ( a : = 0)|^^=0 ( 16 ) 

and for the setup b they will be similar at the right in¬ 
terface as well. However considering the top layer as the 
connecting layer between the right lead and the scatter¬ 
ing region (setup a), at the right interface x = d we must 
have. 


^R{x = d) U, = '^B{x = d)\^^, 

^R{x = d)\B^ =^B{x = d)\g^, 

^B{x = d)\^^=0. ( 17 ) 

Having both scattering matrices Sa and Sb: the gen¬ 
eral condition for bound states as = SaSbclb implies 
that det(l — S^Sb) = 0. Subsequently from Eqs. (2) and 
(7), we find the following characteristic equation for the 
ABS energies, 

del {l - See[S(e) - S*(-£)] - SeeS(£)SeeS*(-e) 

-[1 - SeeS(£)]SfteS(£)[l - SeeS(£)](-£) } = 0, 

( 18 ) 


which can be used to find the Josephson current across 
the junction given by Eq. (1). 



Figure 2. (Color online) The current-phase relation 
/(0)//o(T) in the short junctions which made of undoped BG 
flakes are shown at the zero temperature (T = 0) varying the 
scaled bias voltage Vo /Ao for the first setup in (a) and for the 
second one in (b). The length of the junction is fixed at the 
d/lA = 4. 


III. RESULTS AND DISCUSSION 


Here the numerical results for the Josephson current 
passing through BG flakes between ML regions contacted 
by superconductors are presented. First we assume the 
undoped BG with /i = 0 and study the current phase 
relation (CPR) vary the vertical bias. We denote the 
Fermi wave vector with kp and define l± = hvF/t±_ 
which is a length scale over which the excitations trav¬ 
eling in the two layers of the BG are coupled. Fig¬ 
ure 2 shows the supercurrent variations /(0) scaled by 
/o(T) = 2eA{T)WkF/d for the two different setups when 
the junction length is d/l_\_ = 4 and at the zero tem¬ 
perature T = 0. We see that at the absence of verti¬ 
cal bias Vb ftie setup (a) is in the so-called 7r-state with 
^ sin(7r-h0) while the other is in a 0-state reveal¬ 
ing a CPR of the form /((/))//q ^ sin cj). Very intriguingly 
when the vertical bias is applied the first setup can pass 
a 0 — TT transition when Vb ^ I.SAq as we see in Fig. 2(a). 
However the second setup remains always in 0-state and 
increasing the vertical bias only suppresses the current. 
Moreover we note that the amplitude of current in two 
cases is very different and for the considered parameters, 
setup (b) has almost three orders of magnitude larger 
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Figure 3. (Color online) Critical currents /c//o as function 
of normalized temperature T/Tc are depicted versus different 
values of scaled bias voltage Vo/t±, shifted on the axis, for 
the setup (a). The dimensionless length of the junction is fixed 
at value d/l± =4. The Cusp in the curves represents the tt —0 
transition where the critical current vanish. 


supercurrent. 

It is worth to note that the main difference between two 
setups is the role of interlayer hopping t± in their trans¬ 
port properties even in normal (non-superconducting) 
state. In setup (a) t± plays a significant and constructive 
role by facilitating the pass way through the flake and 
one expect that if we could switch it off the current will 
vanish and transport through the BG regions is not pos¬ 
sible. This could explain why the supercurrent in device 
(a) is much smaller than the device (b) since the elec¬ 
trons must change the layer and enter another transport 
channel. However, in the setup (b) the interlayer hopping 
acts only as a barrier for the movement of quasiparticles 
and it somehow provides a resistance against transmis¬ 
sion through the scattering region. In fact in this case 
only the bottom ML is responsible for the transport and 
subsequently varying none of the parameters like vertical 
bias, the temperature, and doping leads to the qualitative 
change in the behavior of supercurrent and 0 — tt transi¬ 
tions cannot take place. In other words the upper layer of 
the flake only introduces an extra channel for scattering 
off the bottom layer which slightly suppresses the elec¬ 
tron transport through the system. When Vq is absent, 
the device (b) shows large transmission probabilities and 
as a result the CPR is strongly non-sinusoidal while in¬ 
creasing the vertical bias, which opens a gap in the band 
structure, leads to the decline in the supercurrent and 
CPR becomes closer to the sinusoidal behavior. 

In the remainder we will concentrate on the first device 
to understand the underlying physics of supercurrent re¬ 
versal and 0 — TT transitions which can occur by varying 
Vb and /i of the BG flake as well as temperature and the 
junction length when either Vq or ja have a finite value. 


We present the dependence of the critical current Ic/Io 
on the temperature varying the vertical bias scaled by the 
superconducting gap Vb/Ag in Fig. 3. When no vertical 
bias voltage is applied, the junction remains always in tt 
state for all temperatures below the critical temperature 
T < Tc. Nevertheless applying Vb the supercurrent as 
the function of T/Tc shows one cusp indicating a tt to 0 
transition. Moreover the position of the cusps varies with 
Vb and after Vb ^ 1.5Ao the junction will be completely 
in 0 state irrespective of the temperature variations. As 
one can see from Fig. 3 while at small values of vertical 
bias the critical current shows an overall decline, when 
it enters the 0 state for a wide range of temperatures, Ic 
can even increase with T. 

In order to understand above mentioned features es¬ 
pecially the mechanism of supercurrent reversal, we use 
the intuition based on the scattering matrix for the nor¬ 
mal transport governed by Sb which contains the prop¬ 
erties of the BG flake contacted by MG regions. We 
have already discussed that the transmission through 
the junction (a) is very small and subsequently the An¬ 
dreev bound states are formed from single scatterings 
from the junction and the multiple scattering processes 
are strongly suppressed. Upon the electron hole conver¬ 
sion or vice versa at the NS interfaces the quasiparticles 
acquires an energy dependent Andreev phase beside the 
superconducting phase ±0/2. Moreover the excitations 
passing between two interfaces may find an extra phase 
shift 7 sc, corresponding to the phase factors of the trans¬ 
mission coefficients 1^2 and t 2 i defined by t = |t| exp( 7 sc)- 
Therefore one can convince himself that the phase accu¬ 
mulated in the excitations during scattering 7 sc is added 
to the superconducting phase difference which leads to a 
shift in the 0 dependence of the ABS energies. The occur¬ 
rence of TT state for the device (a) when /i = Vb = 0 signals 
that there must be some phase shift of amount tt origi¬ 
nating from 7 sc while in the other setup there is no phase 
shift. The phase 7 sc in the setup (a), in fact, originates 
from the transmission between lower and upper layers of 
the BG flake accompanied with the change of sublattice 
index Ai B 2 . To understand its root let us have a look 
to the Hamiltonian (10) in which the two Dirac Blocks, 
corresponding to the two MLs constructing the bilayer, 
differs with each other. The difference comes from the 
fact that we have AB stacking and the upper layer is 7 r /3 
rotated with respect to the lower. Subsequently the wave- 
functions spinor structures in two layers are not the same 
which leads to the emergence of an extra phase 7 sc when 
the electrons need to pass from lower to the upper layer. 
At this point one may wonder how the two setups become 
different since in both of them the excitations undergoes 
transitions between the layers several times depending 
on the length d. But as one can immediately see from 
Fig. 1 in the second setup since at the end quasiparti¬ 
cles again leave the flake to the lower ML and we have 
the same number of passing from lower to the upper and 
vice versa there will be no net phase accumulated in the 
transmission coefficients. On the other hand in setup (a) 
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Figure 4. (Color online) Critical currents /c//o versus nor¬ 
malized temperature T/Tc and the normalized length of the 
junction djl^ are plotted for the scaled bias values [0.6, 0.9, 
1.2]. The Cusp in the curves represents the tt — 0 transition 
that again moves with the inclusion of topological term. 


there is always one more transition from lower to the up¬ 
per layer and the scattered excitations from the junction 
will have a net phase factor. It must be mentioned that 
7sc in general depends on the angle of incidence of the 
particles, as well as the doping /i and bias Vq. Thus in 
the absence of the bias and doping the two setups are in 
TT and 0 states, respectively, as the overall effect of scat¬ 
tering phases 7sc. It is worth to note that if we had A A 
stacked bilayer flake even in the first setup no phase fac¬ 
tor appears due to the transition between the layers and 
subsequently the junction will be in the 0 state. 

We have already seen that when a vertical bias is ap¬ 
plied, the junction undergoes a tt to 0 transitions which 
is clear from both Figs. 2 and 3. Up to now we have 
considered a junction with fixed length d/l± ^ 4. But if 
we increase the length as shown in Fig. 4 the Josephson 
current changes to a 0 state in a certain length which 
depends on the strength of the vertical bias and the tem¬ 
perature. Moreover at the higher Vb the domain of the 
TT junction versus T/Tc and d/l±_ shrinks and becomes 
smaller. The bias induced transition originates from the 
dependence of the scattering phase 7sc to Vq. In fact 
application of the bias between the layers leads to an 
asymmetry which influences the scattering processes and 
the phase of the different modes transmission coefficients 
such that after a certain value of the vertical bias depend¬ 
ing on the length, the junction becomes of 0 type. It must 
be noted that when the vertical bias voltage goes beyond 
the superconducting gap, the excitations carrying the su¬ 
percurrent becomes evanescent due to the gap formation 
in the spectrum of BG flake. 

The effect of doping is somehow similar to the applica¬ 
tion of Vb and leads to the tt — 0 transitions as indicated 
in Fig. 5. For small values of bias and doping in compar¬ 
ison with superconducting gap the junction remains in 
the TT state while when either /i or Vb proceed well above 



Figure 5. (Color online) Critical currents Ic/Io as function of 
scaled doping of the Bg flake /x/Aq are demonstrated versus 
different values of scaled bias voltage Vo/t±, shifted on the z 
axis, for the setup (a) at the zero temperature T = 0. The 
dimensionless length of the junction is d/l± = 4. The Cusp in 
the curves indicates the tt — 0 transition. 


Aq the supercurrent shows a 0 behavior. In this state 
the critical current becomes larger upon increasing any 
of them. We can relate the transition with fi to the fact 
that by increasing the doping, more modes with differ¬ 
ent incidence angles participate effectively in transport 
and therefore their scattering phases 7sc are washed out. 
This cause the junction to enter a 0 state but due to the 
mechanism different from the bias induced 0 state. 

Finally we should comment on the possible experimen¬ 
tal realization of the device under consideration and the 
predicted results. As it has been mentioned in the in¬ 
troduction, the step-like ML/BL interfaces and bilayer 
flakes connected to monolayer regions have been already 
observed and some transport experiments have been done 
based on them. On the other hand there are a large 
amount of experiments on the superconducting proxim¬ 
ity effect and Josephson junctions based on graphene. So 
it seems that the device we propose here is completely 
reliable in the current experimental setups and the super¬ 
current reversal can be investigated in them by applying 
the vertical bias voltage, changing the doping and even 
the temperature. 


IV. CONCLUSION 

In this paper we employ scattering theory to study the 
supercurrent flowing through a weak link containing a 
biased BG flake embedded between two MG based super¬ 
conducting reservoirs. We investigate two possible con¬ 
figurations for the setup regarding how the BG flake is 
attached to the MG regions in the left and right. The 
reservoirs are connected either to different layers of the 
BG flake or to the same layer of the BG. The supercurrent 
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passing through the undoped BG flake in the first situa¬ 
tion is in the tt state while the other shows a 0 behavior in 
the absence of the bias voltage (Vq = 0 ). Both states are 
robust against varying the length of the junction d as well 
as the temperature T as long as /i = Vb = 0. We argue 
that the underlying mechanism of supercurrent reversal, 
characterized by the cusp in the critical current curves 
in the first setup, is related to the scattering phases ac¬ 
cumulated in excitations upon transmission through the 
junction. We reveal that the biased BG flake in the sec¬ 
ond setup remains always in 0 state and increasing the 
vertical bias only suppresses the current. Nevertheless, 
the biased junction in the first setup undergoes a tt — 0 


transition at zero doping provided by a finite value of Vb 
depending on the length. Increasing the doping causes 
TT — 0 transitions as well so that the junction is in tt state 
as long as // and Vb are small. 
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